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ABSTRACT 

By combining high spatial resolution and wide-field spectroscopy performed, 



> . respectively, with SINFONI and FLAMES at the ESO/VLT we measured the 

ly-v . radial velocities of more than 600 stars in the direction of NGC 6388, a Galactic 

0^ . globular cluster which is suspected to host an intermediate-mass black hole. Ap- 

• . proximately 55% of the observed targets turned out to be cluster members. The 

O ' cluster velocity dispersion has been derived from the radial velocity of individual 

stars: 52 measurements in the innermost 2", and 276 stars located between 18" 

and 600". The velocity dispersion profile shows a central value of ~ 13 km s~^, 

L^ . a fiat behavior out to ~ 60" and a decreasing trend outwards. The comparison 

;_i . with spherical and isotropic models shows that the observed density and velocity 

- - - dispersion profiles are inconsistent with the presence of a central black hole more 

massive than ~ 2OOOM0. These findings are at odds with recent results obtained 

from integrated light spectra, showing a velocity dispersion profile with a steep 

central cusp of 23 — 25 km s~^ at r < 2" and suggesting th e presence of a black 



hole with a mass of ~ 1.7 x IO'^Mq (JLiitzgendorf et al.ll201ll ). We also found some 



evidence of systemic rotation with amplitude A^ot ~ 8 km s ^ in the innermost 
2" (0.13 pc), decreasing to A,ot = 3.2 km s"^ at 18" < r < 160". 

Subject headings: Globular clusters: individual (NGC 6388); stars: evolution - 
black hole physics 
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INTRODUCTION 



The confirmation of tlie existence of intermediate mass (10'^ — 10'^Mq) black holes 
(IMBHs) would have a dramatic impact on a number of open astrophysical problems, rang- 
ing from the formation of supermassive BHs and their co-evolution with galaxies, to the 
origin of ultralum inous X-ray sources i n nearby galaxies, up to the detection of gravita- 
tional waves (e.g.. iGebhardt et al.ll2005l ). However, the evidence gather ed so far in support 



of the existence of IMBHs are in conclusive and controversial (see, e.g., iNoyola et al.ll2010l : 



Anderson fc van der Marelll2010l ). Galactic globular clusters (GCs) are thought to be the 
best places wh ere to search for these e lusive objects. In fact, the extrapolation of the "Magor- 
rian relation" (JMagorrian et al.l 119981 ) down to the IMBH masses naturally leads to the GC 
mass regime. Moreover, numerical simul ations have shown that the co res of GCs are the ideal 
habitat for the formation of IMBHs (e.g. iPortegies Zwart et al.ll2004l ). For these reasons, the 



recent years have seen an increasing number of works dedicated to the search for IMBHs 
in GCs, exploiti ng all observational channels, ranging frorn the detection of X-ray and ra- 
dio emission (see lStrader et al.ll2012l : iKirsten fc Vlemmingsll2012l . and references therein), to 
the detailed study of the cluster st ructure, especially in terms of the shape of the density 



and velocity dispersion profiles (e.g. .IGebhardt et al.l2000t iGerssen et al.ll2002 ilLanzoni et al. 



2011 



hereafter 



20071 : iNovola et al.ll2010l : lAnderson fc van der Marelll2010l : iLiitzgendorf et al.l 
Lll). In spite of such an effort, however, no firm conclusions could be drawn to date. This 
is mainly because of the great difficulties encountered from both the theoretical and the 
observational points of view. In particular, the expected X-ray and radio emission are still 
quite uncertain. A power-law density profile S*(r) ex r" with a slope a ~ —0.3 (significantly 
shallower than the one expected in a post-core collapse system : a ^ —0.7) is expec ted in 



the innermost cluster regions (e.g. iBaumgardt. Makino. &: Hutll2005l: Miocchi 



2007), even 



though this could be a non-univocal signature of a central IMBH (IVesperini fc Trentill2010l ). 
In the presence of an IMBH a steep central cusp is expe cted also in the velocity dispersion 
profile (JBaumgardt. Makino. fc Hutll2005l : iMiocchil 120071 ). However, the precise determina- 
tion of the stellar density and (especially) the velocity dispersion is very challenging in the 
highly crowded central regions of GCs. In particular, measuring the projected velocity dis- 
persion (cTp) in GCs has always been a very difficult task. In fact, GCs are close enough 
for allowing to resolve single stars; hence, the estimate of crp from the line broadening of 
integrated light spectra can be easily falsified by the dominant contribution of a few bright 



^ Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile (under 
proposals 381.D-0329, 073.D-0211 and 073.D-0760). Also based on observations with the NASA/ESA HST 
(Prop. 19835), obtained at the Space Telescope Science Institute, which is operated by AURA, Inc., under 
NASA contract NAS5-26555. 
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objects (e.g., iDubath. Meylan &: Mayor 1119971 ). On the other hand their stellar densities 
(especially in the centre) are so high that the measurement of individual radial velocities is 
very challenging from the observational point of view and it requires very high spatial resolu- 
tion spectroscopy. For the accurate determination of the internal proper motions, both high 



spatial resolution imaging and long enough time baselines are n eeded (e.g. iMcNamara et al. 



2OO3I : lAnderson &: van der Marelll2010l : iMcNamara et al.ll2012l ). Finally, different method 



ologies for measuring the v elocity dispersion could bririg to incompatible results (see the case 
of u Centauri discussed in Ivan der Marel fc Anderson! I2OIOI : iNoyola et al.l 120101 ). 



This paper is devoted to NGC6388. Using high spatial resolution (0.025" x 0.027"/pix) 
data acquired with the HST/ACS High Resolution Channel (HRC), we discovered a shallow 
central cusp {a = — 0.2) in the in ner 1" and, from the comparison with King models including 
a central IMBH (JMiocchil 120071). we suggested that a ~ 6 x IO^Mq BH could be hidden in 
the center of this cluster ( jLanzoni et al.l 120071 ). However, subsequent i nvestigations, both 



in the X-ray and in the radio band s, argued against such a conclusion (JNucita et al.l 12008 



Cseh et al.ll2010l : iBozzo et al.ll201ll ). The question has been recently re-opened by Lll, who 



used integrated light spectra acquired with the FLAMES/ARGUS to derive ap from the 
spectral line broadening: they found an increasing trend toward the cluster centre, up to 
values of 23 — 25 km s~^, and concluded that this is best fitted by assuming the presence 
of an IMBH of (1.7 ± 0.9) x IO^Mq. Here we derive the velocity dispersion profile of NGC 
6388 from the radial velocities of more than 300 individual stars. We find no evidence of 
central cusp in the velocity dispersion and the comparison with theoretical models suggests 
that this is inconsistent with the presence of an IMBH more massive than ~ 2OOOM0. 

The paper is organized as follows. Observations and data analysis are described in 
Section [2l Results are presented in Sect. |3] and include the determination of the cluster 
membership (Sect. 13. ip . the study of its systemic rotation (Sect. 13. 2p and the measurement 
of the velocity dispersion profile (Sect. 13. 3p . Section Hj is devoted to the comparison with 
previous results and theoretical models. The conclusions are presented and discussed in Sect. 
|5l Appendix A is dedicated to the discussion of the impact of unresolved background and 
stellar blends on the determination of radial velocities. 



OBSERVATIONS AND DATA ANALYSIS 



2.1. The SINFONI data set for the cluster central region 

In order to acquire spectra of individual stars i n the center of NGC 6 388 we exploited 
the high spatial resolution capabilities of SINFONI ( lEisenhauer et al.ll2003l ). a near-IR (1.1- 
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2.45 /iin) integral field spectrograph fed by an adaptive optics module and mounted on 
the YEPUN (VLT-UT4) telescope at the ESO Paranal Observatory. By using the 100 
mas plate-scale and the i^-band grating (sampling the 1.95-2.45 fim wavelength range), we 
derived spectra at a resolution R = 4000 for ~ 6 stars located in a 3.2" x 3.2" region 
centered on the cluster gravity center (as quoted in lLanzoni et al.l 120071 ). The observations 
have been performed in service mode between April and June 2008 (ESO proposal ID: 381. D- 
0329(A), PI: Lanzoni), under an average seeing of ~ 0.8" (FWHM). A star of magnitude 
i? ~ 12 located ~ 9" from the cluster center has been used as natural guide star for adaptive 
optics correction, leading to an average Strehl ratiqj of ~ 30%. The high level of crowding 
combined with the small field of view prevented us from estimating the sky background 
through a dithering technique. Hence, for background subtraction purposes the observations 
were performed with a target-sky-sky-target sequence, using an empty sky region located 
~ 2' North-East from the cluster center. A single target frame is the combination of 10 
exposures, each one 20 seconds long. The observations have been repeated until a total 
exposure time of ~ 2.2 hr on target was reached. 

We derived a wavelength calibrated 3D data cube, as well as a reconstructed 2D im- 
age, by using the SINFONI pipeline v.2.3.2 (JModigliani et al.ll2010l ). Following the pipeline 
prescriptions, we first corrected all the target and sky data for darks, fiats, geometrical dis- 
tortions and differential atmospheric refraction. We then subtracted the sky background 
by using the observations of the sky field, we calibrated in wavelength by using a Th-Ar 
reference arc lamp, and we finally combined the corrected target frames to build the 3D 
data cube. In principle this provides a total of 32 x 64 spectra. In practice, in order to 
minimize the contamination of the spectrum from the light of neighbouring stars, we made 
the conservative choice of extracting only the spectrum of the central (and most exposed) 
spaxel corresponding to the star centroid. To this end we perfo rmed the photometric a naly- 
sis of the reconstructed SINFONI image by using ROMAFOT (JBuonanno et al. Ill983l ) and 
thus determ ined the positi o n of t he source centroids. Then, the HST/ACS-HRC catalogue 
discussed in iLanzoni et al.l (120071 ) was roto-translated on the SINFONI one. The resulting 
astrometric solution has an accuracy of less than 0.2 spaxels and provided us with a clear-cut 
identification of 90% of the SINFONI targets. In the remaining cases (corresponding to a 
candidate SINFONI target surrounded by a few close, fainter companions), a careful visual 
inspection of the SINFONI image was necessary to precisely locate the spaxel correspond- 
ing to the brighter star. With IRAF we thus extracted the ID spectrum for a total of 59 



^The Strehl ratio one of the main parameters characterizing the image quahty of data taken with adaptive 
optics systems and it corresponds to the amount of Ught contained in the diffraction-Umited core of the point 
spread function, with respect to the total flux. 
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spaxels (Figure [T]) corresponding to 59 stars in the magnitudes range 13 ^ / ^ 18 (Figure 
|2]). The extracted spectra have been cleaned from the telluric absorption lines through the 
division by a spectrum of a featureless early-type star observed each night with the same 
instrument setup of the science target. The typical signal-to-noise (S/N) ratio per pixel at 
the wavelength 2.18/im ranges from ~ 20 at / = 18, to ~ 80 for the brightest sources. 



2.2. The FLAMES data set for the cluster external regions 



The spectra of individual stars beyond the cl uster core have been acquired by using 
the ESO VLT multi-object spectrograph FLAMES JPasquini et al.ll2002[ ). in the MEDUSA 
UVES-I- GIRAFFE combined mode. This provides 132 fibers to observe an equivalent number 
of targets in one single exposure, over a field of view of 25' in diameter. Under proposal 
381.D-0329(B) (PI: Lanzoni) four pointings of 2320 s each have been performed in June and 
July 2008 with the GIRAFFE grating HR21 (which samples the Ca II triplet spectral range 
at a resolution R = 16200) and the UVES setup Red Arm 580, covering the wavelength 
range 4800A < A < 6800A at a resolution R = 47000. The targets have been selected 
from the comb i ned H ST/ ACS-WFC and ESO/WFI photometric catalog already discussed in 
Lanzoni et al.l (|2007| ) and iDalessandro et al.l ( 120081 ) . considering only isolated objects (with 
no brighter stars within a circle of 1" radius), having l^ < 17 and being located along 
the canonical evolutionary sequences of the color- magnitude diagram (CMD). About 15 
GIRAFFE fibers and 2 UVES fibers in each pointing have been used to sample the sky. A 
total of 422 and 7 stars were observed with GIRAFFE and UVES, respectively. Additional 
MEDUSA data have been retrieved from the ESO Archive. These include the spectra of 
red giant branch stars observed with GIRAFFE gratings HRll and HR13 under program 
073.D-0211 (PI: Carretta), with HR13 under program 073.D-0760 (PI: Catelan), and with 
UVES Red Arm 580 setup in both cases. While some of these stars are in common with the 
sample discussed above, 99 objects are new. The FLAMES data set therefore consists of a 
total of 528 individual star spectra. 

All these data have been reduced with the most updated versions of the UVES and 
GIRAFFE ESO pipelines, including bias subtraction, flat-field correction, wavelength cali- 
bration through a reference Th-Ar lamp, rebinning to a constant pixel-step, extraction of 
the one dimensional spectra and (for the UVES targets only) order merging. The accuracy 
of the wavelength calibration was checked by measuring several emission sky lines in the 
available spectr al ranges and compariri g the observed wavelengths to the rest-frame posi- 
tions quoted by lOsterbrock et al.l (119961 ) . The removal of the sky background is particularly 
crucial for the observations performed with the GIRAFFE setup IIR21, since the third com- 
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ponent (A = 8662A) of the Ca II triplet is heavily affected by O2 emission lines. For targets 
observed with HR11/HR13 and Red Arm 580 setups, instead, the sky level is only a few 
percent of the star one (because of the brightness of the targets), introducing only a neg- 
ligible amount of noise in the stellar spectra. In all cases the sky background has been 
removed by subtracting from each stellar spectrum a master sky spectrum obtained as the 
median of the observed sky spectra. Each individual sky spectrum has been inspected by 
eye in order to check for any possible contamination from close stellar components and a 
few sky spectra showing absorption lines due to the hght of close stars have been discarded. 
Note that the employed pre-reduction procedure takes into account the fiber to fiber relative 
transmission, in order to properly scale the intensity of each spectrum to the different fiber 
response. For the targets observed with HRll and/or HR13 setups, multiple exposures of 
the same stars are available. In these cases we added together the sky-subtracted spectra to 
increase the S/N ratio. The same procedure has been adopted for the UVES targets with 
multiple exposures. The final S/N per pixel is always larger than 50. 



2.3. Stellar radial velocities 

The radial velocities (V^) of the SINFONI targets have been measured from the ^^C^^O 
band-heads (see Figure E]) and, occasionally, a few atomic lines. For the FLAMES targets 
we used several atomic lines, dependi ng on the GIRAFFE s etup. In all cases we adopted 



the Fourier cross-correlation method (JTonry fc Davied Il979l ) as implemented in the fxcor 



IRAF task. The observed spectrum is cross-correlated with a template of known radial 
velocity and a cross-correlation function (i.e., the probability of correlation as a function of 
the pixel shift) is computed. Then, this is fitted by using a Gaussian profile and its peak 
value is derived. Once the spectra are wavelength calibrated, the pixel shift obtained by 
the Gaussian fit is converted into radial velocity. We employed different reference template 
spectra according to the adopted spectral configuration and stellar type. Basically, we used 
suitable synthetic spectra computed in the wavelength range covered by the instrumental 
setup and convolved with a Gaussian profile to reproduce the spectral resolution of the 
used gratings. F or the SINFONI targets we used synthetic spectra computed with the 



code described by lOriglia. Moorwood fc Olival (119931 ) and followi ng the same procedure (i.e. 



model atmospheres, atomic and molecular line list) described in lOriglia fc Rich! ( 120041 ). For 



the F LAMES targets we computed synthetic spectra with the Linux version of the SYNTHE 



code (ISbordonell2005l ). 



Four SINFONI spectra are of low quality. Three of them are among the faintest stars 
in the i^-band and their spectra have a very low S/N ratio (< 5), probably because of a 
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bad response of the corresponding spaxels. The spectrum of the star located at the south 
edge of the field is unusable as a consequence of the adopted jittering procedure. Hence they 
have been excluded from the subsequent analysis. After a careful analysis of the effects of 
unresolved background and blended sources on the determination of the radial velocity (see 
Appendix A), three additional targets have been removed and the final SINFONI sample 
consists of 52 stars located between 0.2" and 1.9" from the cluster center. From the FLAMES 
spectra we measured the radial velocity of all the analyzed targets, corresponding to a total 
of 528 stars located at distances ranging between 11" and 883". 



The velocity errors are computed for each star following the prescriptions by lTonry &: Davies 



( 119791 ). which take into account the accuracy in the fit of the peak of the correlation func- 
tion, the errors due to mismatches between object and template spectra and the effect of 
the anti-symmetric noise component of the correlation function. The typical uncertainties 
are of a few km s~ for the SINFONI targets, 0.4 — 0.5 km s~ for the FLAMES giants, and 
1.5 — 2 km s~ for the FLAMES horizontal branch stars because of the blending between the 
Ca II lines and the Hydrogen Paschen lines. All radial velocities have been reported in the 
same heliocentric reference system, by applying the corresponding correction computed with 
the IRAF task rvcorrect. The presence of several targets in common among the different 
FLAMES datasets allowed to verify possible off-sets in the radial velocity zero point, due to 
the different spectral configurations and wavelength calibrations. We assumed as reference 
sample the one observed with the HR21 setup, since it represents more than 70% of the 
entire dataset. No relevant offset is found when this is compared with the other samples, the 
average radial velocity differences being smaller than 1 km s~^, with a standard deviation 
< 1.5 km s" . This guarantees a good internal homogeneity of our radial velocity zero point 
and confirms the reliability of our error estimates. 



2.4. Stellar metallicities 

With the aim of providing an additional constraint to the cluster membership, we also 
computed the metallicity of our target stars. Unfortunately, a precise measure for the SIN- 
FONI targets is prevented by the low spectral resolution of the available data. We could 
just derive a rough estimate of the target metalhcities from a few Fe, Ca and Na lines, sug- 
gesting that, within a large uncert ainty (^ 0.2 dex), th ey are compatible with the cluster 



mean value ([Fe/H]= —0.44 ±0.01; ICarretta et al.l 120071 ). To measure the metallicity of the 



FLAMES targets different methodologies have been adopted according to the different spec- 
tral configurations used. The atmospheric parameters (effective temperature and gravity) 
have been estimated for each target from the comparison between the star position in the 



CMD and theoretical isochrones by iPietrinferni et al.l (|2006[ ). For the stars observed with 
gratings HRll and HR13 and with UVES, we performed a direct estimate of the metaUicity 
from the observed equivalent widths (EWs) of the Fe I lines. Given the target atmospherical 
parameters and the instrumental spectral resolution, we selected o nly the Fe I lines predicte d 
to be unblended. EWs were measured by means of DAOSPEC (IStetson fc Pancindl2008l ). 
Details about the choice of t he atomic data and the computation of model atmospheres and 
chemical abundances are in iMucciarelli et al.l (120121 ). In the case of targets observed with 
the HR21 grating, with the exception of very cold stars and blue horizontal branch stars 
(for which the Ca II lines are heavily affected by strong TiO molecular bands and Hydrogen 
Paschen lines, respectively), t he metaUicity has bee n derived from the Ca II lines by follow- 
ing the approach described in ICarrera et al.l (120071 ). The metaUicity has been measured for 
a total of 508 stars. 



RESULTS 



3.1. Cluster membership 



The radial velocities measured from the FLAMES spectra are plotted as a function of 
the distance from the cluster center (r) in the top- left panel of Figure HI Most of the observed 
stars clearly clump around K- ~ 80 km s~ , while the remaining ones show a large spread at 
lower, mostly negative values. Hence, a selection in Vr has been used as driving criterion to 
distinguish cluster members from field objects. 

As apparent from the bottom panels of Fig. HJ the separation between cluster mem- 
bers and field stars is less clear in terms of metaUicity. This is indeed expected, since the 
f ield population is constituted by disk giants with [Fe/H] ranging between zero and — LO 
( JRobin et al.ll2003l ). Hence, to determine the average metaUicity of NGC 6388 we selected 
(280) stars with 50 < V^ < 130 km s~^; in addition, based on the star distribution in the 
([Fe/H], Vr) plane, we excluded the most evident outliers by applying a preliminary cut in 
metaUicity: —0.7 <[Fe/H]< —0.1 (see the dotted box in the bottom-left panel of Fig. H]). 
The average val ue turns out to be ([F e/H]) = —0.40 with a dispersion of 0.08 dex, in good 
agreement with ICarretta et al.l (120071 ). We then conservatively considered as cluster mem- 
bers only those stars having a metaUicity value within 3a from the meano this yields to 



•^Note that by estimating the atmospheric parameters from a theoretical isochrone appropriate for NGC 
6388 (see Sect. I2.4|) we imphcitly ass umed that a l l star s belong to the cluster. While this is not necessarily 



true, the Besancon Galactic Model (|Robin et al.l l2003l ) shows that the expected contamination from field 
stars in the color, magnitude and radial velocity ranges covered by our targets is quite low (^ 3.5%). 
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276 stars (54% of the total) with —0.64 <[Fe/H]< —0.16, which are marked as sohd circles 
in the bottom- and top-left panels of Fig. HI Their metallicity distribution is compared to 
that of the entire population (508 objects) in the bottom-right panel of the figure. No signif- 
icant changes are found for reasonable modifications of the assumed cuts in radial velocity 
and metallicity. Moreover, no significant impact is expected on the following analysis by 
the potential inclusion of a few interlopers. The resulting 276 FLAMES cluster members 
are highlighted as solid circles in Figures |5] and [6] which show, re spectively, the clus t er ma p 



and CMD obtained from the photometric catalog presented in iDalessandro et al.l (120081 ). 
The radial velocities of SINFONI targets range between 56 and 104 km s^^. By also taking 
into account their very central radial position (r < 2"), we consider all of them as cluster 
members. 

To determine the systemic velocity (Kys) of NGC 6388 we adopted even more con- 
servative limits: out of the 276 cluster members, only (240) stars with r < 350" and 
60 < Vr < 105 km s^ were selected (see the solid box in the top-left panel of Fig. 
H]). The mean value, computed with a Maximum-Likelihood method, turns out to be 
V^ys = 82.0 ± 0.5 km s~^. Different, but still reasonable, assumptions about the limits 
in cluster-centric distance and/or radial velocity cuts do not produce any significant change 
in this result, which is also in agreement with the s ystemic velocity quoted in th e literature 



for NGC 6388 (see, e.g., Harris 1996 - 2010 version; ICasetti-Dinescu et al.ll2010f ). The coor- 
dinates, magnitudes and radial velocities (referred to l^yg) measured for all (52-1-276) cluster 
members are listed in Tabled] (where we also give the metallicities of FLAMES targets). 
In the rest of the paper, V^. will indicate radial velocities referred to the cluster systemic 
velocity: Vr = Vr — Vsys. 



3.2. Systemic rotation 

Figure [7] shows l^ as a function of the distance from the cluster center for the 52 
SINFONI targets. A bimodality is apparent in this distribution and suggests the presence 
of systemic rota t ion. T o investigate this possibility, we used the method fully described in 



Bellazzini et al.l (120121 . and references therein). We considered a line passing through the 



cluster center with position angle PA varying between 0° (North direction) and 90° (East 



Moreover, adopting an isochrone appropriate for the field contaminating population ([Fe/H]=0, distance =8 
kpc), the resulting metallicities turn out to be only a few 0.1 dex larger. In any case, we stress that the 
driving selection criterion for cluster membership is the radial velocity, while metallicity is used just to refine 
the selection, even at the cost of rejecting a few cluster members. Most importantly, we stress that any 
metallicity selection is unable to introduce kinematical biases in the analysis. 
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direction), by steps of 15°. For each value of PA, such a hue sphts the SINFONI sample 
in two. The difference A{Vr) between the mean radial velocity of the two sub-samples was 
computed and it is plotted as function of PA in the upper panel of Figure El Its coherent 
sinusoidal behavior is a signature of rotation and the parameters of the best-fit sine function 
provide us with the position angle of the rotation axis (PAq = 11°) and the amplitude of the 
ro tation curve (Arnt = 8 -5 km s~ ; for the exact meaning of this parameter see the discussion 



m 



Bellazzini et al.ll2012l ). The corresponding rotation curve is shown in the lower- left panel 
of Fig. [HI The Kolmogorov-Smirnov probability that the distributions of Vr for the two 
sub-samples on each side of the rotation axis are drawn from the same parent population is 
~ 4% (lower-right panel), indicating a not particularly high statistical significance, possibly 
due to the limited sample. We verified that taking into account or ignoring such a rotation 
has only very weak effects on the computation of the cluster velocity dispersion, which is 
dominant. While a rotation or shearing signature within ~ 3" from the centre is claimed 
also by Lll, it seems to be mainly caused by the shot noise produced by three bright stars in 
the field of view (see their Figure 7). In any case, our best-fit position angle of the rotation 
axis seems not to coincide with their. 

We performed a similar analysis also for the FLAMES sample, considering various dis- 
tance intervals from the cluster centre. The strongest signal of systemic rotation is found 
for stars with 18" < r < 160" (160 in total). The corresponding parameters are PAq = 200° 
and Afot = 3.2 km s~^, and the Kolmogorov-Smirnov probability for a common origin of 
the Vr distributions for the two sub-samples is of only 0.6% (see Figure [H]). No evidence of 
significant rotation is instead detected in the outermost cluster regions, where the second 
moment of the velocity distribution decreases and it becomes increasing difficult to measure 
independently the or dered and the r a ndorn motions. These results do not confirm those 



recently published by iBellazzini et al.l (120121 ). who suggest that some rotation is present be- 
tween ~ 160" and 400" from the center, based, however, on a small sample of less than 40 
stars (in the same radial range we count, instead, more than 100 stars). 

We finally emphasize that the results presented in Lll are all in terms of the second 
moment of the velocity distribution (I^ms = \/Kot + '^pj with V^rot indicating the rotational 
velocity). Hence, a direct comparison with their work (Sect. 14. ip requires to analyze the 
same quantity. Also for this reason, exactly as in Lll, we do not apply any correction for 
rotation and we deal with the second moment of the velocity distribution. However, because 
the rotation is small, in the following we continue to discuss in terms of velocity dispersion 
and we use symbol crp instead of I^ms all over the paper. 
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3.3. Velocity dispersion 

The radial distribution of Vr for all the selected cluster members is shown in Figure 
[TOl To compute the projected velocity dispersion profile crp(r), the surveyed area has been 
divided in a set of concentric annuli, chosen as a compromise between a good radial sampling 
and a statistically significant number ( ^ 50) of stars. In each radial bin, o"p has been 
computed from the dispersion of the values of K measured fo r all the stars i n the annulus, 



by following the Maximum Likelihood method described in I Walker et al.l (l2006l . see also 
Martin et al. 2007; Sollima et al. 2009). An iterative 3cr clipping algorithm was applied in 
each bin (the stars thus excluded from the computation are marked as grey circles in Fig. 
[TU1) . The error estimate is performed by following Pryor & Meylan (1993). The resulting 
velocity dispersion profile is shown in Figure [11] and listed in Table [2l Given the number 
of stars in the SINFONI data set (52 objects in total) we computed o"p(r) by considering 
both one single central bin (solid square in the figure) and two separate annuli (26 stars 
at r < 1.2", 26 stars beyond; see the empty squares). For the external data set we tried 
different sets of radial bins (grey region in Fig. [TTl) and the values obtained by adopting the 
four annuli listed in Table [2] are shown as black circles. These results show that the velocity 
dispersion of NGC 6388 has a central (r ~ 1") value of ~ 13 km s~ , stays approximately 
flat out to r ~ 60", and then decreases to ~ 7 km s~^ at 200" < r < 600". If no a clipping 
algorithm is applied, the velocity dispersion profile remains almost unchanged, with the only 
exception of the outermost data-point that rises to ~ 8.9 ± 0.8 km s~^. 



4. DISCUSSION 

4.1. Comparison with previous studies 

The inner part (r < 10") of the velocity dispersion profile derived by Lll from the line 
broadening of integrated-light spectra (triangles in Fig. [TT]) is clearly incompatible with the 
one obtained in this study from the radial velocities of individual stars (squares). While the 
exact reason for such a disagreement is not completely clear, a detailed comparison between 
our radial velocity measurements and Lll radial velocity map (which is available in electronic 
form and shown in their Figure 7) suggests that the shot noise corrections applied by these 
authors may have been insufficient. 

Figure \T2\ approximately reproduces the central 2" x 2" region of Lll's Figure 7, with 
the values of Vr that we measured for each individual star in the same field of view marked in 
black. In the second annulus, for instance, we measure radial velocities as low as 73, 75, 77 km s~ 
in regions where they quote V^ > 98 km s~^ (see the red spaxels on the right-hand side). 
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We also find Vr = 88, 99 km s~^ where they measure K < 74 km s^^ (cyan region on the 
upper-left side). This might be because those spaxels in Lll are still highly contaminated, 
respectively, by the very bright stars with Vr > 100 km s~^ and Vr ^ 62 km s~^ (in both 
studies) that produce the large red spot on the right-hand side of the figure and the dark-cyan 
spot at the top, despite the applied shot noise correction (see the white asterisks flagging 
a few spaxels in these red and dark cyan regions). More in general, we note that, even in 
the case of close stars, our individual radial velocity values can show significant differences, 
thus indicating that the measurements are independent each other. Instead, the Vr color 
codes of Lll are quite uniform in several spaxels around the brightest stars (i.e., the de- 
gree of correlation, or cross-talking, is high). Even though Lll do not use the values of Vr 
to measure the cluster velocity dispersion (which is estimated from the line broadening of 
the combined spectra), this analysis suggests that the spectra of the spaxels close to the 
brightest stars are still significantly contaminated and provide an artificial broadening of the 
combined lines. On the other hand, our estimate is directly derived from the dispersion of 
the radial velocities measured for each individual star in the field of view and the statistics 
(52 stars) is quite good. Hence we do not see what kind of mistake we could have committed 
to artificially decrease the value of crp(r). We finally note that the center adopted by Lll 



I 



is slightly offset from ours (see Fig. [T2l) Ji However, even if we adopt their center and their 
radial bins, we still find results that are totally inconsistent with a steep central cusp: we 
measure ap = 14.5 km s~ in the central bin (r < 0.9", 17 stars), ap = 12 A km s~ in the 
second annulus (0.9" < r < 1.9"; 27 stars), and ap = 13.1 km s~^ if a single bin at r < 1.3" 
(44 stars) is adopted. 

The deriv ed central ve l ocity dispersion also di sagrees with the value o f 18.9 ±0.8 km s~ 
determined by llUingworthl (119761 ) and reported by lPryor fc MeylanI (119931 ). We beheve that, 
even in this case, the difference comes from the fact that such a value has been obtained from 
the line broadening of integrated-light spectra, which is prone to shot noise bi as. However , we 
also note that v alue of 18.9 km s~^ refers to the inner 6" — 12" from the centre IlllingworthI (see 
Sect. II in 119761 ) . while the effective velocity dispersion quoted by Lll has been computed 
for r < 40". T his cast doubts up on the claimed agreement between the value quoted by 
Lll and that of IlllingworthI (119761 ). In addition, the latter author attempted to correct the 
observed velocity dispersion, thus to derive the central value; the result is quoted in his 
Table 5: 19.6 ± 5.7 km s~ . Hence, within such a large uncertainty, both our measurement 



''Since there is a typo in the coordinates of the cluster centre quoted in equations (2) and (3) of Lll, we 
have used their Figure 6 to approximately locate their centre on the HST/ACS-HRC image. As shown in 
the left-hand panel of Fig. [1^ our centre is located 0.55" North and 0.40" East from theirs (differently from 
what is stated in Lll). 
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and that of Lll a re formally in agreement with the central velocity dispersion obtained by 
lUingworthl JiQTgI ). 



4.2. Comparison with models 



In this section theoretical models are used to reproduce the observations and derive 
some constraints on the presence of a central IMBH in NGC 6388. To this end, we follow 
two different and complementary approaches. First, starting from a family of self-consistent 
models admitting a central IMBH, we select the one that best reproduces the observed 
density and velocity dispersion profiles. This provides us with the corresponding structural 
and kinematical parameters in physical units, including the IMBH mass. Second, we use the 
observed density profile and include a variable central point mass, to solve the spherical Jeans 
equation (this is what is done also in Lll). We therefore obtain the corresponding family 
of velocity dispersion profiles, which are then compared to the observations to constrain 
the BH mass. The first approach is more realistic in terms of the stellar mass-to-light 
ratio [M/L), since it takes into account various populations of stars with different masses 
and different radial distributions (as it is indeed expected and observed in mass segregated 
GCs), while in most applications of the Jeans approach a constant M/L profile is assumed 
for the computation of the velocity dispersion (but see, e.g., Williams et al. 2009 and 
Liitzgendorf et al.ll2012l . for Jeans models with varying M/L). Moreover the self-consistent 
modeling has the advantage that the models stem from a distribution function which is known 
a priori. Hence, the model consistency (i.e., the non-negativity of the distribution function in 
the phase-space) is under control and guaranteed by construction. However, the validity of 
the results is limited to the case in which the assumed models are the correct representation 
of the true structure and dynamics of the cluster. The second approach is m ore general. 



but the resulting models could be non-physical (e.g. iBinney &: Tremaind 119871 ). In all cases 
we assume spherical symmetry and velocity isotropy. T hese are reaso nable assumptions for 
NGC 6388, which does not show significant el lipticity (JHarris Ill996l ) and is thought to be 



quite dynamically evolved fJFerraro et al.l 120121 . see also Lll 



Self- consistent King and Wilson models — In order to reproduce the observed surface density 
and velocity dispersion profiles we used self-consistent, m ulti-mass King (1966) and Wilson 
(1975) models admitting the presence of a central IMBH ( Miocchil 120071 ) ii We assumed the 



■5 As detailed in iMiocchil |2007h . the IMBH is included following the iBahcall fc Wolj |l976l) phase-space 
distribution function. These models can be generated and freely downloaded from the Cosmic-Lab web site 
at the address: http: //www. cosmic-lab. eu/Cosmic-Lab/Products .html 
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same stellar mass spectrum as in iLanzoni et al.l ( 120071 ). consisting of six mass bins, O.lMg 



'0 



wide and ranging between O.SMq and 0.9Mq, plus a seventh bin containing a population 
of 1.2Mq white dwarfs, assumed to be the remnants of 4-8 Mq stars. For the comparison 
with the observations we used only the density and velocity dispersion profiles of the [0.8, 
0.9] Mq mass group, corresponding to turnoff and giant sta rs (the ones eff ectively used in 



the observations). By following the procedure described in iMiocchil ( l2010l ). we performed 
a parametric fit to the observed density profile: this provided us with the model structural 
parameters (gravitational potential, characteristic radius and ratio between the BH mass 
and the cluster total mass). Once the projected density is fitted, the shape of the velocity 
dispersion profile is univocally determined and we therefore applied only a vertical shift 
(corresponding to a velocity scale factor) to adjust it to the observations. As shown in 
Figure [121 both the King and the Wilson models properly reproduce the observed profiles, l 



Once the distance of NGC 6388 is specified (we adopted d = 13.2 kpc from lDalessandro et al. 



20081 ). the assumed velocity scale factor univocally determines the value of the cluster total 
mass (Mtot)- In turn, given the BH-to-cluster mass ratio set by the density best-fit, this 
provides us with the mass of the IMBH. The resulting values are listed in Table [3l with M^^ 
indicating the BH mass obtained from the self-consistent modeling. This approach suggests 
an IMBH of ~ 2OOOM0. 

Jeans models — Starting from the best-fit King and Wilson density profiles (upper panels of 
Fig. [T3|) . we solved the Jeans equation by admitting the presence of an IMBH with a mass 
(Mgjj) ranging between zero and thirty times the best-fit BH mass of the corresponding 
self-consistent model (i.e., M^h/^bh = 0^ --'SO, with Mfn = 2147M0 and 2125M0 for the 
King and Wilson cases, respectively; see Table |3]). ResultCl are shown in Figure [H] and 
indicate that our observed velocity dispersion profile is consistent with both the absence of 
a central IMBH, and the presence of a black hole with a mass up to ~ 2000Mq. Instead we 



^Just for illustrative purposes, we show in Fig. [10] the ±3cr velocity dispersion profiles corresponding to 
the best-fit King and Wilson models (dotted and dashed curves, respectively). As apparent, only a few stars 
used in the computation of ap{r) as described in Sect. 13.31 are located below or above these lines, and the 
velocity dispersion profile that one would obtain by excluding them from the analysis (consistently with the 
assumption that the stellar system behaves as a perfect King or Wilson model) is still in agreement with the 
one quoted above. 

^For fixed BH mass the same velocity dispersion profile is expected in both approaches. However, the 
comparison between (Jp(r) shown in Fig. [T3l and the curve corresponding to ^-^^ehM^bh = 1 in Fig. [T4| 
reveals a slightly different shape. This is because the profile shown in Fig. [T4| corresponds to the velocity 
dispersion produced by all cluster stars, instead of that produced by the [0.8, 0.9] Mq mass group only 
(which is shown in Fig. I13p . In fact, the adopted Jeans equation approach is unable to distinguish among 
different mass groups (or, in other words, it assumes a constant M/L ratio). 
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find that in order to match the Lll profile an IMBH of ~ 6 x lO^M© should be assumed, in 
rough agreement with the value quoted by these authors. 

The results of our analysis are consistent with both no IMBH and a black hole of no more 
than 2OOOM0 in the centre of NGC 6388. In particular, its presence is mainly constrained by 
the cusp of the density profile, while its mass is set by the observed velocity dispersionjfl In 
fact, the steep central increase expected in o"p(r) is not visible in our observations, because, 
for a 2OOOM0 BH, it is predicted to occur at a very small distance from the centre (r ^ 0.5" = 
0.03 pc), where we measured only 9 stars. Note that such a small cusp radius is in agreement 
with the standard computation of the BH sphere of infiuence (reH ~ GMbh/ct"^, with o 
being the unperturbed velocity dispersion): in fact, assuming a ~ 14 km s~ , one gets 
tbh ~ 0.7". We finally conclude by noticing that the discussed family of self-consistent 
models are incompatible with the velocity dispersion profile derived by Lll (triangles in Fig. 
[T3!) . In fact, in these models the velocity dispersion cusp due to the presence of an IMBH 
always appears within the innermost 3" — 5" from the centre and the main effect of increasing 
the BH mass is to make the cusp steeper, while there is no mean of moving it outward, at 
10" or more, where Lll already observed a significant increase of o"p(r). 



5. SUMMARY AND CONCLUSIONS 

By using the high spatial resolution spectrograph SINFONI, and the multi-object spec- 
trograph FLAMES/MEDUSA installed at the ESO VLT, we measured the radial veloci- 
ties and metallicities of ~ 600 stars in the direction of the Galactic globular cluster NGC 
6388. Based on the measurements of VJ. and [Fe/H] obtained for the FLAMES sample, we 
determined the systemic radial velocity and mean metallicity of NGC 6388: Kys = 82.0 ± 
0.5 km s~^ and ([Fe/H]) = — 0.40±0.08 (which are in good agreement with the values quoted 
in the Literature). We also set the criteria for cluster membership: 60 < K < 105 km s~ 
and —0.64 <[Fe/H]< —0.16. This provided us with 52 members located at 0.2 < r < 1.9" 
and 276 between 18" and 600" from the SINFONI and the FLAMES samples, respectively. 

The radial velocities of cluster members were used to study the systemic rotation and 



8 



Indeed, the difference in th e derived value of Mb h with respect to our previous estimate (5700Mo) is 



mostly due to the fact that in iLanzoni et al.l (J2007h we used a photometric esti mate of the clust er total 



mass under the common assumption of a stellar mass-to- light ratio M/Lv — 3 (e.g. lDiorgovskilll993l ). Here, 



instead, the measured velocity dispersion provides us with a dynamical estimate of the cluster total mass, 
which turns out to be a facto r of ^ 2.2 smaller. This is in agreement with recent results suggesting that 



M/Lv -- 1.3 in GCs (see, e.g. ISollima et al.ll2012l) 
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velocity dispersion of NGC 6388. We found evidence of systemic rotation with amplitude 
^rot ~ 8 km s~^ in the innermost 2" (0.13 pc), decreasing to Aj-ot = 3.2 km s~^ at 18" < r < 
160". In the outermost regions no rotation signature is detected, most probably because the 
sample is not rich enough and the (largely dominant) velocity dispersion hides the rotation 
signal. The cluster velocity dispersion profile, both in the very center (0.2" < r < 1.9") and 
in the external regions (out to ~ 600") has been derived from the dispersion of individual 
star radial velocities. Results indicate a central value of ~ 13 — 14 km s~^, a flat behavior 
out to ~ 60" and a decreasing trend outwards. No evidence of a central cusp in crp(r) is 
detected in our data. 

These results desagree with those published by Lll, who find a steep cusp at r < 10" 
and a central velocity dispersion of ~ 23 — 25 km s~^. The reason for such a discrepancy 
is not completely clear, but it is likely related to the fact that Lll measure up from the 
line broadening of integrated light spectra. As it is widely discuss ed in the literature (e.g., 
Dubath. Meylan fc Mayor 1119971 : lAnderson fc van der Marell 120101 ). this procedure is prone 
to strong biases in the case of Galactic GCs. In fact the acquired spectra can be dominated 
by the contribution of a few bright stars, instead of effectively sampling the underling stel- 
lar distribution. Certainly, using high-resolution IFU spectroscopy, applying procedures to 
correct for the shot noise induced by giants, and measuring the line broadening from the 
combined spectra of several spaxels in the same radial annulus (as it is done in Lll) may 
alleviate the problem. However, (i) the fact that it is hard to imagine methodological errors 
able to artificially lower the value of crp(r) if estimated from individual radial velocities (at 
least if the statistics is good, as it is the case in our study which comprises 52 stars), (ii) 
the apparently high degree of color correlation among several spaxels close to the brightest 
stars in the Vr map of Lll, and (Hi) the disagreement (in terms of both velocity dispersion 
and radial velocity map; see Figs. [11] and [12]) between the two methodologies suggest that 
measuring cxp from the integrated light spectra can be fallacious. Indeed, similar discrep- 
ancies between the "integrated light" and the "statistical" approaches are found also in the 
case of u Centauri, where the latter method is appli ed to proper motion measurements (see 
Anderson fc van der Marell boiol : iNovola et al.lboiOl ). 



The derived velocity dispersion profile is consistent with the presence of an IMBH at 
most as massive as ~ 2OOOM0. Given the measured central velocity dispersion (13-14 
km s~ ), the derived IMBH mass is in agreeme nt with the expectat ions of the Mbh — cr 
relation within its quoted uncertainties (see e.g., JGiiltekin et al.ll2009l ). More sophisticated 
models are needed to put stronger constraints. However, we believe that to obtain a firm 
answer it is first necessary to precisely understand why integrated light spectra and the 
measurement of individual star radial velocities (or proper motions) yield to different val- 
ues of the velocity dispersion. It is also necessary to collect more data, providing the full 
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velocity profi le (rotation, velocity d ispersion, and also the higher order terms; see, e.g., the 



discussion in iGebhardt et al.ll2005l ). To this end, high spatial resolution photometry and 



spectroscopy are key observational tools, and new soph isticated methods for the analys of 



integral field unit data cubes (e.g.. iKamann et al.ll2013l ) are very promising. Meanwhile the 



James Webb Space Telescope and the planned ground-based extremely large telescopes will 
become available, for a significant progresses in this field of the research it is certainly worth 
to push even further the usage of the HST and the current generation of spectrometers at 
8- 10m class telescopes, with integral field capabihties and with adaptive optics. 

We warmly thank the referee for the careful reading of the manuscript and the useful 
suggestions that helped improving the quality of the paper. This research is part of the 
project COSMIC-LAB funded by the European Research Council (under contract ERC- 
2010-AdG-267675). 



A. IMPACT OF UNRESOLVED BACKGROUND AND BLENDS ON THE 

DETERMINATION OF K 

We estimated that the contamination from unresolved background stars (that, in prin- 
ciple, could artificially lower the derived velocity dispersion) is totally negligible. In fact, the 
dominant contribution to the background is due to unresolved / ^ 19 stars which have effec- 
tive temperatures Tes ~ 5500 K. On the other hand, for Tes > 4500 K the depth of the CO 
band-heads rapidly decreases for increasing temperature (see Figure [TB]) . and at T^s = 5500 
K it is totally negligible with respect to the CO depth at 5000 K (the temperature of our 
faintest target). We therefore conclude that the unresolved background does not affect the 
determination of Vr for our targets. We also stress that, even if, we took into account only 
the brightest targets (45 stars with / < 16.5), the resulting central velocity dispersion would 
remain essentially unchanged: up = 13.8 ± 1.5 km s~^. 

In the case of non-isolated targets, we estimated the effect of blends on the computation 
of Vr as follows (note that more complete methods have been proposed in the literature. 



but t hey are far too much sophisticated for the quality of our data; see, e.g., iKamann et al. 



20131 ). By convolving each star with its point spread function, we estimated the fraction 
of light from neighbouring objects that contaminates the spaxel of the candidate SINFONI 
target. We found that the contribution of stars detected in the HST image but not corre- 
sponding to a SINFONI centroid was always negligible (< 5%). This is not surprising, since 
these stars are fainter by at least 1-1.5 magnitudes than the SINFONI targets. Instead, we 
found seven couples of SINFONI targets where the spaxel corresponding to the fainter star 
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(and used to measure its K; see Sect. I2.3P was contaminated by 20-25% from the light of the 
neighbouring (brighter) star. For each of these couples, we simulated the spectrum of the 
fainter object blended with that of the neighbouring star, by taking into account the percent- 
age of contamination and the differences in flux, temperature and radial velocity between 
the two "companions". We then measured the radial velocity from this blended spectrum 
and we compared it to the input one (i.e., the value previously measured from the observed 
spectrum of the fainter target). In four cases, the resulting value of Vr was in agreement 
with the input one within the uncertainties. Discrepancies of about —5 km s~^ have been 
found only for the remaining three cases for which the contaminating (bright) star had Vr 
smaller than the contaminated (fainter) one. Such a behaviour is due to the asymmetry of 
the CO band-heads: when the radial velocity of the contaminating star is larger than that of 
the fainter target, the offset is hidden within the CO band shape, otherwise it is detectable 
as a shift in Vr towards lower values. We therefore conclude (in agreement with the findings 



of iKamann et al.l 120131 ) that blending is not a major issue for the measured radial velocities 
but, to be conservative and to avoid any possible concern about the velocity measurements, 
we decided to remove the three most affected targets from the sample. Note that including 
or excluding these stars (or even all the seven targets with larger contamination) has no 
impact on the derived velocity dispersion. 
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id 


RA 


DEC 


B 


V 


/ 


Vr 


et/r 


Fe/H] 


data set 


10000131 


264.0715464 


-44.7349731 


- 


14.94 


12.82 


23.5 


2.3 


- 


1 


10000216 


264.0720428 


-44.7353238 


- 


15.12 


12.77 


-18.1 


2.0 


- 


1 


10000027 


264.0712349 


-44.7352065 


- 


15.53 


13.80 


13.5 


2.0 


- 


1 


10000325 


264.0724705 


-44.7351295 


- 


15.53 


13.78 


-18.7 


2.7 


- 


1 


30828 


264.1096197 


-44.7367217 


16.90 


14.96 


— 


-4.8 


0.3 


-0.49 


2 


7000555 


264.0853010 


-44.7688900 


16.90 


15.27 


13.40 


-13.8 


0.2 


-0.48 


2 


7001285 


264.0194000 


-44.7440720 


17.34 


15.98 


14.37 


1.1 


0.1 


-0.39 


2 


7001303 


264.1128920 


-44.7231290 


17.74 


16.32 


14.61 


10.0 


0.2 


-0.41 


2 



Table 1: Identification number (id), coordinates (in degree), magnitudes in the B, V, and 
/ bands, radial velocity and its uncertainty in km s~'^ (V^ and e^^, respectively), [Fe/H] 
abundance, and flag for the data set (fl=l meaning SINFONI, fl=2 indicating FLAMES) for 
all the stars selected as cluster members and used to compute the rotational velocity and 
velocity dispersion of NGC 6388. No values of metallicity are given for the SINFONI targets 
because of their large uncertainty. The complete version of the table is available online. 
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Table 2: Velocity dispersion profile of NGC 6388. The three first columns give the internal, 
external and mean radii (in arcseconds) of each considered radial bin (r^ is computed as 
the average distance from the centre of all the stars belonging to the bin), A^^ is the number 
of star in the bin, ap and Ca-p are the velocity dispersion and its rms error (in km s~^), 
respectively. 
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Model Wo r, Mtpt Mfn 

King 11 6.4 1.2 x 10^ 2147 

Wilson 11 6.8 1.2 x 10^ 2125 

Table 3: Parameters of the King and Wilson models that best fit the observed surface density 
and velocity dispersion profiles (see Figure [T3|) . Wq is a dimensionless parameter proportional 
to the gravitational potential at the cluster center, r^ is the core radius in arcseconds (defined 
as the radius at which the projected stellar density drops to half its central value), Mtot and 
Mglj (in units of Mq) are, respectively, the total cluster mass and the mass of the central 
IMBH estimated from the self-consistent modelling approach. 
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Fig. 1.— Reconstructed SINFONI image. The field of view is 3.2" x 3.2", Nortli is up, East 
is left. Red crosses mark the spaxels from which we extracted the spectra, corresponding to 
stellar centroids measured in the HST/ACS-HRC image. Green circles mark the four stars for 
which the radial velocity could not be measured and the three targets conservatively rejected 
because significantly blended with a close brighter star (see Sect. 12.31 and App endix A). The 



large blue cross flags the adopted position of the cluster gravity centre (from iLanzoni et al. 
2007h . 
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Fig. 2. — HST/ACS-HRC color-magnitude diagram of the central region of NGC 6388. The 
59 stars for which we extracted a SINFONI spectrum are highlighted as large circles: solid 
circles correspond to the 52 objects for which we could reliably measure the radial velocity, 
empty circles flag the remaining seven stars. 
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Fig. 3. — Comparison between one spectrum acquired with SINFONI and the appropriate 
template. Three CO band-heads used to compute the stellar radial velocities are flagged. 
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Fig. 4. — Radial velocity and metallicity distributions (top and bottom panels, respectively) 
of the FLAMES targets used to determine the cluster systemic velocity, mean metallicity and 
membership. The mean metallicity of NGC 6388 (see label in Panel d) has been estimated 
from the (280) stars enclosed in the dotted box (50 < Vr < 130 km s~ and —0.7 <[Fe/H]< 
—0.1) drawn in Panel c. Only the (276) stars having a metallicity value within 3a from the 
mean have been considered as cluster members (solid circles in Panel c). Their metallicity 
distribution is compared to that of the entire population in Panel d (shaded and empty 
histograms, respectively). The radial velocities of the 276 cluster members (solid circles) are 
shown as a function of the distance from the cluster centre in Panel a. The corresponding 
number distribution is plotted as shaded histogram and compared to that of the entire 
FLAMES sample (empty histogram) in Panel b. The box drawn in Panel a encloses the 
sub-sample of members used to conservatively compute the cluster systemic radial velocity, 
which is labelled in Panel b. 
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Fig. 5. — Map of stars brighter than V = 18, derived from the combined HST/ACS-WFC and 
ESO/WFI photometric catalog discussed in IPalessandro et all J2008I ). The 276 FLAMES 
targets selected as cluster members are highlighted as large solid circles. The concentric 
annuli used to compute the cluster velocity dispersion profile in the external regions (see 
Table |5]) are also marked. 
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Fig. 6.— {V, B-V) color-m agnitude diagram of NG C 6388 obtained from HST/ACS-WFC 
and ESO-WFI observations JPalessandro et al.lboosh . The 276 FLAMES targets selected as 
cluster members are highlighted as large solid circles. 
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Fig. 7. — Left panel: Distribution of K- — V^ys for the 52 SINFONI targets. The histogram 
in the right-hand panel shows a bimodahty indicating the presence of systemic rotation. 
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Fig. 8. — Upper panel: Difference between the mean radial velocities of the SINFONI targets 
located on each side of the cluster with respect to a line passing through the center and having 
position angle PA (measured from North, PA=0°, to East, PA=90°), as a function of PA. 
The solid line is the sine function that best fits the observed pattern. The corresponding 
values of the rotation amplitude (Ay-ot) and the position angle of the rotation axis (PAq) are 
labelled. Lower left panel: Distribution of Vr — K)ys as a function of the distance X' from 
the center projected onto the axis (dashed line) perpendicular to the best-fit rotation axis. 
Lower right panel: Cumulative radial distribution of stars having X'(PAo) > (solid line) 
and X'(PAo) < (dotted line). Pks quotes the Kolmogorov-Smirnov probability that the 
two distributions are drawn from the same parent population. 
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Fig. 9.— As in Fig. El but for the 160 FLAMES targets located between 18" and 160" from 
the cluster centre. 
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Fig. 10. — Radial velocity (referred to the cluster systemic velocity) as a function of the 
distance from the center, for the selected cluster members: 52 SINFONI targets are shown as 
crosses, 276 FLAMES stars are plotted as solid circles. The radial bins (see Table |2]) adopted 
to determine the external portion of the velocity dispersion profile from the FLAMES data 
are marked as rectangles. Their vertical size corresponds to the 3cr range of the local velocity 
dispersion crp(r); stars not included in these boxes (and not used in the computation of ap) 
are marked as grey circles. The dotted and dashed curves mark the ±3cr velocity dispersion 
profiles of the best-fit King and Wilson models, respectively (see Sect. 14. 2p . 
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Fig. 11. — Projected velocity dispersion profile of NGC 6388 as determined from the radial 
velocity of individual stars (squares and circles). The solid square marks the value of o"p 
obtained from the entire SINFONI sample (52 stars at r < 1.9"), the two empty squares 
correspond to the values derived within and beyond r = 1.2" (26 and 26 stars in the two bins, 
respectively). Results from the FLAMES dataset are shown as black circles, with the grey 
region indicating the dispersion obtained for different choices of the radial bins. Triangles 
correspond to the results obtained by Lll from integrated light spectroscopy and are shown 
for comparison. 
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Fig. 12. — Left panel: Approximate reproduction of LH's Figure 7, showing in colors their 
radial velocity map for the innermost 2" x 2" (North is right, East is up). The values of Vr 
increase from ~ 30 km s~^ (dark blue), to more than 100 km s~^ (red), as indicated in the 
color-bar on the right-hand side. White asterisks mark the spaxels that Lll excluded from 
the analysis to correct for the shot noise effect due to the presence of very bright stars. The 
values of the radial velocities that we measure for each individual star detected in the same 
field of view are marked in black at the star position. The two circles (of radius 0.9" and 
1.9") correspond to the two innermost radial bins considered in Lll. Right panel: /-band 
HST/ACS-HRC image of the same cluster region. The orientation of the map and the two 
circles are as in the left-hand panel, while the dashed box indicates the field of view of our 
SINFONI observations. Red crosses mark the stars for which we measured the radial velocity 
Vr, which is labelled in black in the left panel of the fig ure. The blue cross m arks the position 
of the gravity centre used in the present work (from iLanzoni et al.l 120071 1 . while the green 
circle indicates the cluster centre adopted in Lll. 
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Fig. 13. — Comparison between the observations and the best-fit self-consistent King and 
Wilson models with central IMBH (left- and right-hand panels, respec t ively) . The observed 
density profile (solid circles in the upper panels) is from lLanzoni et al.l (120071 ). The observed 
velocity dispersion profile (symbols in the the lower panels) is the same as in Fig. [TTJ Solid 
lines show the best-fit model results for the mass group corresponding to 0.8 — 0.9Mq stars. 
Model parameters are given in Table |3l 
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Fig. 14. — Comparison of the observed velocity dispersion profile (the same as in Fig. [TTj) 
with two families of Jeans models. Left panel: solid lines correspond to Jeans models calcu- 
lated from the King density profile shown in the upper- left panel of Fig. [13] and by assuming 
different black hole masses: from bottom to top, M^^/M^^i = 0) 0-5, 1, 2, 3, 5, 10, 30, with 
Mfjj = 2147M0. Right panel: the same as in the left panel, but for Jeans models calculated 
from the Wilson density profile shown in the upper-right panel of Fig. [131 In this case, 
M^H = 2125M0. 
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Fig. 15. — Synthetic spectra in the main CO band-head region plotted for five different 
values of the effective temperature (from bottom to top: Tcs = 3500,4000,4500,5000,5500 
K). 



